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1. Introduction 

Lattice QCD simulations are currently being performed with dynamical degenerate u- and 
d- quarks with a mass close to their physical value as well as the strange quark, using a number 
of different discretization schemes with the most common being Wilson-improved, staggered and 
chiral fermions. Furthermore, simulations at several lattice spacings and volumes are becoming 
available, enabling a comprehensive study of lattice artifacts. The masses of low-lying hadrons 
have been computed and extrapolated to the continuum limit using large enough lattice sizes to 
ensure that volume effects are small [1, 2]. These calculations show agreement with experiment 
and therefore pave the way for evaluating other phenomenologically interesting quantities beyond 
these masses. 

Several collaborations, using dynamical quarks with pion mass down to about 300 MeV, have 
calculated the pion electromagnetic (EM) form factor [3] , which is obtained from the matrix ele- 
ment (7t + (p')\Ji 1 \n + (p)) = + p'n)F K (q 2 ), where q 2 = (p' — p) 2 = —Q 2 . Based on vector dom- 
inance, lattice data are fitted to the form F 7[ (Q 2 ) = (l + (r 2 )<2 2 /6) 1 to extract the mean squared 
radius, which is shown in Fig. 1. As can be seen, there is an increase in the value of (r 2 ) at small 
pion mass, m n . An accurate extraction of (r ) benefits from evaluating the form factor at small 
values of Q 2 accomplished by using twisted boundary conditions (b.c). In a recent calculation, 
ETMC combined twisted b.c. and the so called 'one-end' trick to incorporate the all-to-all propa- 
gator and improve statistics. Using simulations with two degenerate light quarks (Nf = 2) at two 
lattice spacings and two volumes [4] the assessment of cut-off and volume effects was carried out. 
Lattice results on F n obtained with pion masses in the range of 300 MeV to 500 MeV, are extrap- 
olated to the physical point using NNLO chiral perturbation theory (PT). The resulting form factor 
is shown in Fig. 2 [4] and it is in agreement with experiment. 
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Figure 1: The pion mean square radius as a function of Figure 2: F K extrapolated to the physical pion 
ml obtained using simulations with N f = 2 twisted mass mass (red band) using NNLO chiral PT corn- 



quarks. 



pared to experiment (blue band). 



As simulations with quark masses close to the physical value become available, the study of 
resonances and decays of unstable particles becomes an important issue. The p -meson width has 
been studied by several groups [5]. Considering a n + TC~ system in the / = 1 -channel, the P-wave 
scattering phase shift 8\ \ (k) in infinite volume is related via Liischer's relation to the energy shift in 
a finite box. Using Nf = 2 twisted mass fermions (TMF) and considering the center of mass frame 
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and two moving frames one extracts the phase shift at different values of the energy, shown in 

Spnrt k? 



Fig. 3. From the effective range formula tan5n(&) 



6^ E CM (Ml-El M ) 



, where k 



/4- 



one determines Mr and the coupling g p7l7[ and then extracts the width using F p = -g^^, where 



Icr = \JM\j\ — m\. The results on the width as a function of m 2 are shown in Fig. 4 [6]. 
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Figure 3: The p -meson phase shift at m K 
308 MeV for a lattice of L = 2.8 fm. 



Figure 4: The p-meson width for Nf = 2 twisted 
mass fermions as a function of ml.. 



Having reproduced the low-lying hadron spectrum [1, 2, 7, 8], the masses of excited states can 
be studied using e.g. variational methods [9]. Furthermore, one can go beyond masses and consider 
form factors (FFs) and generalized partem distributions (GPDs) that probe hadron structure. The 
characterization of nucleon structure, in particular, is considered a milestone in hadronic physics 
and many experiments have been carried out to measure nucleon FFs and structure functions. Ex- 
periments on nucleon FFs started in the 50s. A new generation of experiments using polarized 
beams and targets are yielding high precision data spanning larger Q 2 ranges. Therefore, nucleon 
FFs serve as a further benchmark for lattice QCD. FFs provide ideal probes of the charge and mag- 
netization distributions of the hadron as well as a determination of its shape in analogy to similar 
studies in e.g. deuteron and other nuclei. 
Non-relativistically the form factor can be related 
to the density distribution via 

F(q 2 ) = fd 3 xe-%* < V|P(*)|V >• 
In Fig. 5 we show the intrinsic charge density contours 
of a spin-zero nucleus showing deformation revealed 
through measurements of transition densities using 
electron scattering. 



Figure 5: Tomographic view of the zero-spin de- 
formed nucleus l54 Gd derived from its rotational 
bands using electron scattering. 

In what follows we will review the status of lattice QCD calculations on baryon form factors 
and nucleon generalized parton distributions. 
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2. Nucleon Generalized Form Factors - Definitions 



In this section, we briefly define the quantities for which results are presented. High energy 
scattering can be formulated in terms of light-cone correlation functions. Considering one-particle 
states \p') and \p), GPDs are defined by [10, 1 1]: 

1 fdX 



F r (x^,q 2 ) = -j—e^( P '\y{-Xn/2)r&e 



A/2 

ig J dan A(na) 

¥&n/2)\p), 



where P = (p' + p)/2, I; = — n ■ q/2, x is the momentum fraction, and n is a light-cone vector with 
Pn = \. 

There are three different types of operators, depending on the choice of T. 

Considering nucleon states these are "Handbag" diagram 
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Expansion of the light cone operator leads to a tower of local twist-2 operators ff^ 1 '"^', related to 
moments. The diagonal proton matrix elements (P\ 6y(x) \P) , measured in deep inelastic scattering, 
are connected to the parton distributions q{x), Aq(x), 8q{x). The twist-2 operators are defined by 
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where q = q^ + q^,Aq = q^ — q^,8q = qj +q±, and the curly brackets represent a symmetrization 
over indices and subtraction of traces. The off-diagonal matrix elements extracted from deep virtual 
Compton scattering can be written in terms of generalized form factors (GFFs), which contain both 
form factors and parton distributions: 



(N(p',s')\ffp^\N(p,s)} = u N (p',s') 



£ [A n+u (q 2 )y^+B n+hi (q 



i'=0, even 



2niN 



q^.-.q* 



P^> . . .F M +mod(n,2)C„ +1 . (^ 2 )— qfrq* . . 

m N 



u N (p,s) (2.1) 



and similarly for ff^ (in terms of A n i(q 2 ), B n i{q 2 )) and ffn^ow (in terms of A T nb Bf tj , Cj ; and D^). 

We list the following special cases: 

• n = 1: Ordinary nucleon form factors: 

A w (q 2 )=F l (q 2 ) = j\xH(x,^q 2 ), B w (q 2 ) = F 2 (q 2 ) = J dxE(x,£ ,q 2 ) 
A w {q 2 )=G A (q 2 ) = j\xH{x^,q 2 ), B [0 (q 2 ) = G p (q 2 ) = J dx£(x,£ ,q Z ) , 
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where in the case of the EM current, = \j/(x)y ll yf(x), the nucleon matrix element is written in the 
form %(///) 7 M Fi {q 2 ) + '-^F 2 (q 2 ) u N {p,s). The Dirac F x and Pauli F 2 FFs are related to the 

electric and magnetic Sachs FFs via the relations: Gsiq 2 ) = F\{q 2 ) — jS^n Fi (q ) and Gm(^ 2 ) = 
F\{q 2 ) +F2(q 2 ). For the axial vector current A^ = y(x)YhY5y Y{x) the nucleon matrix element is 
of the form u N (p',s') Y»Y5G A (q 2 ) + j^G p (q 2 ) \u N (p,s). 

• A„o(0), A„o(0), Aj (0) are moments of parton distributions, e.g. {x) q = A2o(0) and (x)a^ = A2o(0) 
are the spin independent and helicity distributions. Knowing these quantities one can evaluate the 
quark spin, J q = \ [A2o(0) + 620(0)] = jAZ^ + L q and investigate the fraction of the spin carried by 
quarks and its contribution to the total spin via the nucleon spin sum rule, ^ = ^AL q + L q + J g , as 
well as the momentum fraction carried by gluons via the momentum sum rule: (x) g = 1 — A2o(0). 

3. Lattice evaluation 

In order to extract the matrix elements connected to GFFs we need to evaluate three-point 
correlators and compute the renormalization of the operators involved. Despite recent progress on 
the evaluation of disconnected loops, most lattice calculations of GFFs do not take into account 
disconnected contributions. Therefore, in what follows, we consider iso-vector operators for which 
such contributions are zero in the isospin limit. For one-derivative operators, mixing with lower di- 
mension operators is avoided by symmetrizing over the Lorentz indices and making them traceless. 
The study of cut-off and finite volume effects in a systematic way has just begun for baryon GFFs. 
These are more difficult to assess since chiral expansions that describe such dependencies are not 
as developed as in the light meson case. The presence of more uncertainties in the chiral expansion 
combined with the larger statistical noise, which for the nucleon two-point function increases like 
ip^e ^ e (m N -?>m K /2) J ^/yy^ ma ^Q the extrapolation of these quantities to the physical point much 
more demanding. In this review we will focus on: i) Nucleon form factors and lower moments 
using dynamical simulations with pion mass m n ~ 500 MeV and spatial lattice length L ~ 2 fm 
and ii) the N-A system in order to determine the complete set of coupling constants needed in 
chiral expansions. Other topics relevant to hadron structure, such as the strange nucleon FFs, hy- 
peron, Roper and nucleon negative parity FFs, distribution amplitudes and transverse momentum 
dependent parton distributions can be found in Ref. [12] and in contributions to this volume. 

• Three-point functions: For the extraction of matrix elements of local operators we need the 
evaluation of two-point and three-point functions defined by o v 

G(p,t) = £ e- iS fPY\ a (J a (x f ,t f )Jp(0)) " 

x f — — — — (aw) 




G^(F,p',q,t) = £ e^e-vre T pa (J a (x> f ,t f )^ \x,t)J p (x i ,t i )) 

Only the displayed connected diagram is evaluated, which, for most current applications, is done 
by using sequential inversion "through the sink" fixing the sink-source separation tf — 1[, final 
momentum p' and T- projection matrices. Smearing techniques are crucial for improving ground 
state dominance in three-point correlators and thus keep tf — tj as short as possible. We stress that 
it is important to ensure that the time separation tf — 1\ used is sufficiently large by performing the 
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calculation at a bigger time separation and checking the consistency of the results. The generalized 
eigenvalue method can further improve identification of the ground state [13] and enlarge the upper 
range of accessible <2 2 -values [14]. 

• Renormalization constants: Most collaborations use non-perturbative renormalization. Us- 
ing a momentum dependent source [15] one evaluates 



S"(p) = y^e-^-yHu(x)u(y)),G(p) = a ^ £ (u{x)u{z) f {zJW)d(y)) (3.1) 

x,y x,y,z,z' 



with the amputated vertex functions given by T(p) = (S l '(p))~ G(p) (S (p))~ and Jf determines 

the operator, e.g. ^(1,1') = 5 z ,z'V^ D v ' would correspond to the local vector current. The Z- 
factors can be determined in the RI'-MOM scheme by imposing the following conditions 
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' q 12 



1 , (3.2) 



to extract the renormalization factors Z q and Z&. These conditions are imposed in the massless the- 
ory and therefore a chiral extrapolation is needed. The mass-dependence is very weak for the vector 
and axial vector operators. This is demonstrated in Fig. 6 for the case of the one-derivative vector 
the axial- vector operators using Nf = 2 TMF [16]. Results on Zy and Z4 as a function of (ap) 2 
are shown in Fig. 7 where plateaux are improved after subtracting Gia 2 ) -terms perturbatively [17]. 
Using the RI'-MOM scheme but with a momentum independent source, the RBC-UKQCD Col- 
laborations made a comparison between perturbative and non-perturbative determination of the 
renormalization constants and found that results for (x) u -d [18] with perturbative renormalization 
are lower bringing them in agreement with LHPC's results, which used perturbative renormaliza- 
tion [19]. It is therefore important to compute the renormalization constants non-perturbatively. 
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Figure 6: Mass dependence of the renormalization 
constant for vector and axial vector one-derivative 
operators for Nf = 2 TMF [16]. 
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Figure 7: Zy and Za with perturbative subtraction 
of ^(a 2 )-terms for AT/ = 2 TMF [17]. 



• Cut-off effects: The nucleon axial charge the isovector momentum fraction < x > u -d= 
A2o(0) and helicity fraction (x)A»-Arf = ^2o(0) are calculated directly at Q 2 = requiring no fits. 
We can examine their dependence on the lattice spacing by obtaining these quantities at a given 
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value of the pseudoscalar mass in units of ro. In Fig. 8 we show results at three lattice spacings 
using Nf = 2 TMF. As can be seen, G (a 2 ) -terms are small and, allowing a linear dependence, yields 
consistent results to those obtained with a constant fit. This is also true for the nucleon isovector 
anomalous moment K v , Dirac and Pauli radii r\ and r\ that require fits to the EM form factors. 
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Figure 8: Left: gA, (x) u -d> an d (x)au-m\ Right: the nucleon anomalous magnetic moment fc,., the Dirac ra- 
dius r\ and Pauli radius r? times K v as a function of (afro) 2 . The red line is the result of fitting to a constant; 
the blue one is a linear fit. The results are obtained using Nf = 2 TMF [20]. 

We therefore conclude that cut-off effects are small for a < 0.1 fm for ^(a)-improved actions and 
that one can use continuum chiral perturbation theory to extrapolate to the physical limit. 

• Finite volume corrections: In Fig. 9 we compare results on gA, (x) u -a and (x)A U -Ad com- 
puted on different lattice sizes as a function of m 2 
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Figure 9: g A , 3(x) u - d and |(x) A „-Arf using TMF [20], 
Clover [21], hybrid [19] and DWF [18]. 



The observations are: i) Accurate lattice data 
by LHPC using domain wall valence quarks on 
staggered sea (hybrid) for m % ~ 350 MeV with 
Lm n =4.5 and Lm^ = 6.2 show no volume ef- 
fects; ii) TMF results for m n ~ 300 MeV with 
Lm % = 3.3 and Lm % = 4.3 are consistent; hi) 
Results for {x) u -d using Clover fermions from 
QCDSF for m K ~ 270 MeV with Lm n = 3.4 and 
Lm % = 4.2 are consistent, whereas gA differs by 
about a standard deviation; iv) RBC-UKQCD 
results with domain wall fermions (DWF) with 
Lm % = 3.9 and Lm % = 5.7 show no volume ef- 
fects for (x) u - d and (x) A „_ Arf [18]. 



Within the current statistical uncertainties this comparison, therefore shows that volume effects on 
< x > u -d an d < x > A U -Ad are negligible for Lm K ~ 3.3. For gA volume effects of about one standard 
deviation are seen for Lm^~3.4 and pion mass below 300 MeV and therefore larger lattices keeping 
Lm n ~ 4 may be needed. 

Comparing the isovector nucleon EM FFs as a function of Q 2 at m K ~ 300 MeV we also find 
consistent results for Lm n = 3.3 and Lm n = 4.3. A similar behavior is also observed for the nucleon 
axial FF Ga(Q 2 ) whereas the induced pseudoscalar FF G P {Q 2 ), which has a pion pole behavior, 
may suffer from larger finite volume corrections at low <2 2 -values [22]. 
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4. Results on nucleon form factors 

Having examined cut-off and volume effects we compare results from various collaborations 
computed with dynamical fermions for a ~ 0.1 fm 1 and Lm n ~ 3.3. 

• Nucleon axial charge: The axial charge is well known experimentally. Since it is deter- 
mined at Q 2 = there is no ambiguity associated with fitting the <2 2 -dependence of the FF. In 
Fig. 10 we show recent lattice results using TMF, DWF and a hybrid action of DWF on a staggered 
sea, all of which are renormalized non-perturbatively. As can be seen, there is a nice agreement 
among different lattice discretizations and no significant dependence on the quark mass down to 
about m % = 260 MeV. 

1 1 1 1 

* 

• TMF a=0, estimate L-correction 
— TMF at 3 0-values 

I I I I 

"'".00 .05 .10 .15 .20 "'".00 .05 .10 .15 .20 

m, 2 (GeV 2 ) mj 2 (GeV 2 ) 

Figure 10: Left: Lattice data on g A using N f = 2 TMF [20] (a = 0.089 fm: filled red circles for L = 2.1 fm 
and filled blue squares for L = 2.8 fm; a = 0.070 fm: filled green triangles for L = 2.2 fm; a = 0.056 fm: 
purple star for L = 2.7 fm and open yellow square for L = 1.8 fm), Nf = 2+1 DWF [23] (crosses for 
a = 0.114 fm and L = 2.7 fm) and Nf = 2 + 1 using DWF and staggered sea [19] (a = 0.124 fm: open 
orange circles for L = 2.5 fm and open cyan triangle for L = 3.5 fm). The physical point is shown by the 
asterisk. Right: Volume corrected TMF results extrapolated to the continuum limit together with the fit using 
HB^PT (blue band). The band bounded by the lines is the resulting fit to the TMF data shown on the left. 

To illustrate the size of lattice artifacts and obtain a value of gA at the physical point, we use 
TMF results [22]. The volume corrected [24] data are extrapolated to a = using three lattice 
spacings, namely a = 0.089 fm, 0.070 fm and 0.056 fm, at two values of the pseudoscalar mass, by 
fitting to a constant. For intermediate masses we use data at the two coarser lattices. The continuum 
volume-corrected results are shown in Fig. 10. Chiral extrapolation using one-loop heavy baryon 
chiral perturbation theory (HB^PT) in the small scale expansion (SSE) [25] with three fit param- 
eters produces a value of gA = 1.12(8) at the physical point, which is lower than the experimental 
value by about a standard deviation. The large error is due to the strong correlation between the 
A axial charge g^A and the counter-term involved in the fit. Therefore, a lattice determination of 
gAA will allow a more controlled chiral extrapolation. Fitting the raw lattice data produces the band 
shown by the dotted lines. Therefore one observes that, although the continuum volume-corrected 
results are closer to experiment, the largest uncertainty is due to the chiral extrapolation and at the 
physical point the values obtained using the raw and continuum volume-corrected lattice data are 
consistent. 

'We note that results by the LHPC using a hybrid action have a = 0.124 fm. 
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• Nucleon form factors: Recent lattice results on the EM isovector and axial FFs are shown in 
Fig. 11. We observe a nice agreement among lattice results, in particular for Ge(Q 2 ) and Ga(Q 2 )- 
However, both Ge{Q 2 ) and Ga(Q 2 ) decrease with Q 2 less rapidly than experiment. We note that 
a good description of the Q 2 — dependence for both Ge{Q 2 ) and Gm{Q 2 ) is provided by a dipole 
form using the lattice-computed p— meson mass. Lattice results on G P (Q 2 ) using TMF and those 
obtained using the hybrid action on a larger volume, are not consistent, in particular at small Q 2 
where G P (Q 2 ) increases rapidly due to the pion-pole behavior. From the observed quark mass 
dependence of G P (Q 2 ) [22] the 50 MeV difference in the pion mass may not be sufficient to fully 
account for this discrepancy, which may indicate volume effects. 




Figure 11: Left: Isovector electric and magnetic nucleon FFs at n% ~ 300 MeV using TMF [20] (filled blue 
squares) DWF [26] (crosses), hybrid [19] (open orange circles) and Clover fermions [27] (yellow stars). 
Experimental data are shown with the filled green circles accompanied with Kelly's parametrization shown 
with the dashed line. Right: Axial nucleon FFs. The solid line is a dipole fit to experimental data for Ga{Q 2 ) 
combined with pion pole dominance to get the solid curve shown for G p (Q 2 ). 

Using HB^PT to one-loop, with A degrees of freedom and iso-vector N-A coupling included in 
LO [28] we perform a fit to F\ {m K , Q 2 ) and Fj{m % , Q 2 ) with five parameters, namely the iso-vector 
magnetic moment at the chiral limit kJ?, the isovector and axial N to A coupling constants and two 
counterterms. As can be seen, the chiral extrapolation increases the value of F\ and F2 at low Q 2 , 
bringing it into qualitative agreement with experiment. Application of twisted b.c. with a study 
of the associated volume corrections [29] will be very useful in enabling us to obtain these FFs 
at lower <2 2 -values, permitting a better chiral extrapolation. Using the parameters extracted from 
the fits to F\ and F2 we obtain the chiral dependence of the isovector nucleon anomalous magnetic 
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moment and the Dirac and Pauli radii shown in Fig. 12. 
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Figure 12: Left: Open squares show the chirally extrapolated results at the physical point. The dashed line 
is Kelly's parametrization of the experimental data. Right: The solid lines show the prediction of HB^PT 
using the parameters obtained from fitting F\(m K ,Q 2 ) and F2(m K ,Q 2 ). The notation is the same as that of 
Fig. 10. For r\ results using Clover fermions [21] are also shown with the cyan cross-in-square symbols. 



5. Results on nucleon moments 

In this section we show results on the nucleon matrix element of the one-derivative operators 
u Y{ii D v } u — dyip D v } d and uy^y^ D v } " — ^7s7{p D v } d in the MS scheme at a scale ju = 2 GeV. 
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Figure 13: Recent results on the isovector A 20 = (x)u-d and A20 = (x)au-Ai1- 

In Fig. 13 we compare recent results from ETMC [30], RBC-UKQCD [18], QCDSF [21] and 
LHPC [19] on the spin-independent and helicity quark distributions. All collaborations except 
LHPC use non-perturbatively computed renormalization constants. As already mentioned, The 
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ETMC has, in addition, subtracted 0{a 2 ) terms perturbatively to reduce lattice artifacts [16]. There 
is a spread in the values of the lattice results. It was noted that taking a renormalization free ratio 
leads to a better agreement among lattice data with Lni K > 4 [31]. In particular this brought the 
LHPC data in agreement with those from ETMC and QCDSF. 

In HB^PT [32] the expressions for the -dependence of A 2 o and A20 are given by: 
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Using A 2 = 1 GeV and the TMF results we obtain the curves shown in Fig. 13, which yield a 
value higher than experiment for both observables. The very recent result by QCDSF [21] at m % ~ 
170 MeV remains higher than experiment and highlights the need to understand such deviations. 



6. A/-A system 

• Ny* — > A form factors: There is an extensive experimental program to study the N to A 
EM transition and in particular to obtain accurate results on the sub-dominant quadrupole FFs 
G* E2 {q 2 ) and G* C2 (q 2 ) that probe deformation. The experimental results, shown in Fig. 14, are 
compatible with the blue band obtained assuming deformation in the N-A and incompatible with 
the red band that includes no deformation. These FFs can be computed within lattice QCD and 
since no disconnected contributions are involved they provide yet another benchmark for lattice 
methods. The matrix element for N to A EM transition is written in terms of three Sachs FFs as: 



( A{p\s') \f\N(p,s)) = i^u^ G {p\s') 



i\ (Q 2 ) K Ml + G *E2(Q 2 ) K E2 + G *C2 



(Q 2 )K ( 



C2 



where £/ = yf {niNin^/ E^{p')En{p)) 1 ^ 2 is a kinematical factor. 




Hybrid a=0.124 fm m n =353 MeV 
DVF a=0.114 fm m^=330 MeV 
DVF a= 0.084 fm m =297 MeV 



'4^ 



1.0 

Q 2 



1.5 

(GeV 2 ) 



u N (p,s) 




0.4 0.6 

Q 2 (GeV 2 ) 



Figure 14: to A EM transition: Left: The transverse-longitudinal response function Gu vs cm. angle 
between p and 7* (from MAMI and Bates) [33]; The to A magnetic dipole FF (middle) and the ratio of 
Coulomb quadrupole to magnetic dipole FF (right) for the hybrid action and Nf = 2 + 1 DWF 

The extraction of the sub-dominant quadrupole FFs is enabled by constructing optimized 
sources that isolate G* E2 and G* C2 [34] . In Fig. 14 we show results using a hybrid action of DWF on 
staggered sea as well as Nf = 2 + 1 DWF, provided by RBC-UKQCD. For the dominant dipole FF 
G* m , like for the nucleon FFs, we observe a weaker <2 2 -dependence as compared to experiment, 
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that again highlights the need of studying these FFs using simulations with m % < 300 MeV. Lattice 
results, shown in Fig. 14 for the ratio of the Coulomb quadrupole to the magnetic dipole FF are 
non-zero. This is also true for the electric quadrupole FF, pointing to a deformation of the iV-A. 

• Weak N to A transition: In contrast to the EM transition, the weak N to A is not well 
studied experimentally. Therefore a lattice determination of the axial vector N to A FFs would 
provide important input for phenomenology and in particular for chiral perturbation expansions. 
The weak N to A matrix element (A(p',y)|Au|iV(p,j)) = u\(p' ,s')0'x ll u N {p,s) with 



m N 



r 



Ci(Q 2 )_ p , v 



nn 



C A ((1 2 ) 



mill 



where Cj (Q 2 ) is the equivalent of the nucleon FF Ga(Q 2 ) and C^(Q 2 ) of G P (Q 2 ) showing a pion 
pole behavior [34]. In Fig. 15 we show results on the dominant FFs Cj and Cg obtained using the 
hybrid action at m % ~ 350 MeV and with Nf = 2 + 1 DWF at m n ~ 330 MeV and m % = 300 MeV. 



1 ■£ 




q quenched Wilson, m =410 MeV. 
q hybrid, m = 353 MeV 

dynamical Wilson, m = 384 MeV 



Q 2 (GeV 2 ) 



Q 2 (GeV 2 ) 



l GeV 



Figure 15: Left and middle: Axial to A FFs C$ and Cg respectively. The squares are for the hybrid action 
at nij ~ 350 MeV, the filled red circles for DWF at m K ~ 330 MeV and the open blue circles for DWF at 
m K = 300 MeV. Right: A electric quadrupole FF for quenched, Nf = 2 Wilson and Nf = 2 + 1 hybrid action. 

7. A electromagnetic form factors and structure 

Experimentally the A FFs are very difficult to measure due to the fact that the A decays strongly. 
Only its magnetic moment is measured experimentally albeit with a large error. Therefore lattice 
calculations can complement experiment by providing these FFs. 




Figure 16: Quark transverse charge 
densities in the A ++ (left) and D. 
(right) for the 3/2-spin projection 
along the x-axis. Darker colors de- 
note smaller values and the charge 
of the particle is taken into account. 
A dotted circle of radius 0.5 fm is 
included for comparison. 



The matrix element (A(p', s') \f (0) \A(p, s)) = w A ,« {p\s')0 a ^u^ (p,s) can be written as 



q a qP 



F^(Q 2 )g ap +F 4 *(Q 2 ) 



iO^qv 
2m\ 
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with e.g. the quadrupole FF given by: Gei = (F* — tF 2 *) — i(l + t) (F 3 * — tF 4 *), where X = 
Q 2 /(4m 2 J. Optimized sources are constructed to isolate the quadrupole FF Gei, which probes 
deformation. The transverse charge density of a A polarized along the x-axis can be defined in 
the infinite momentum frame. Using Gei we can predict the 'shape' of A. The result is shown in 
Fig. 16 and for spin projection 3/2 it is elongated along the spin axis. The Q. shows a similar but 
smaller deformation [35]. The weak A FFs can be computed in an analogous manner [36]. 

8. Conclusions 

The nucleon EM form factors provide a benchmark for lattice QCD beyond hadron masses. 
Most collaborations obtain results for the isovector FFs up to about Q 2 = 2 GeV 2 . Systematic 
studies of lattice artifacts on GFFs are now under way and recent data reveal that cut-off effects 
are negligible for a ~ 0.1 fm, whereas finite volume corrections, although difficult to evaluate, are 
within the current statistical errors of ~ (2-3)% for Lm n ~ 3.3. A possible exception is G p at 
low <2 2 -values. We find that, in general, lattice results using different discretization schemes are 
consistent but they show a milder <2 2 -dependence as compared to experiment. As illustrated in 
the case of the nucleon axial charge, the biggest uncertainty in comparing with experiment is the 
chiral extrapolation. Therefore a lattice determination of a number of couplings used as input in 
chiral extrapolations will enable global fits to e.g. the N-A system that can help extrapolation to the 
physical point. Interesting questions such as the 'shape' of a hadron can be addressed using input 
from lattice form factors as demonstrated for the A and Q. Moments of GPDs yield more detailed 
information on both longitudinal and transverse distributions and a tomography of hadrons can be 
obtained by studying these quantities. We therefore, conclude that, overall, there is good progress 
in baryon structure calculations and that we now are in an exciting era, having simulations close 
enough to the physical point, in order to probe interesting dynamics in hadronic systems. 
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